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Abstract 

In the last ten years, a number of "Conventional Fluctuation Theorems" have been 
derived for systems with deterministic or stochastic dynamics, in a transient or 
in a non-equilibrium stationary state. These theorems gave explicit expressions 
for the ratio of the probability to find the system with a certain value of entropy 
(or heat) production to that of finding the opposite value. A similar theorem for 
the fluctuations of the work done on a system has recently been demonstrated 
experimentally for a simple system in a transient state, consisting of a Brownian 
particle in water, confined by a moving harmonic potential. In this paper we show 
that because of the interaction between the stochastic motion of the particle in water 
and its deterministic motion in the potential, very different new heat theorems are 
found than in the conventional case. One of the consequences of these new heat 
Fluctuation Theorems is that the ratio of the probability for the Brownian particle 
to absorb heat from rather than supply heat to the water is much larger than in the 
Conventional Fluctuation Theorems. This could be of relevance for micro/nano- 
technology. 
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1 Introduction 

As is well known the First and Second Law of Thermodynamics only involve 
averages of the physical quantities of macroscopic systems but say nothing 
about their fluctuations. In particular the Second Law for Irreversible Pro- 
cesses states that the average ^ heat Q internally produced in an irreversible 
process has to be positive. In the last ten years a number of fluctuation the- 
orems have been derived for the fluctuations of thermodynamic properties in 
non-equilibrium stationary [1,2], as well as transient states [3,4], which consti- 
tute refinements of the laws of thermodynamics in so far that they take into 



Here, as well as in the rest of the paper, a bar means an average. 
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account a property of their fluctuations which goes far beyond the statement 
that Q > 0. In this respect a generahzation of thermodynamics including 
fluctuations has been in progress. 

It should be pointed out that these new fluctuation properties are valid for 
large fluctuations around a non-equilibrium stationary state possibly far from 
equilibrium. As such they differ from fluctuations dealt with in the context of 
Irreversible Thermodynamics, such as the Fluctuation Dissipation Theorem 
and the Onsager relations (as treated e.g. in the classical book by De Groot 
and Mazur [5]), which refer to small fluctuations around equilibrium. 

The Fluctuation Theorems (FT) will be divided into two classes: Conven- 
tional (CFT) and New (NFT) Fluctuation Theorems, where in the latter case 
the word "theorem" is premature. We will not conflne ourselves here to heat 
fluctuations alone but also consider work and energy fluctuations. In the lit- 
erature until now the overwhelming number of papers has dealt exclusively 
with heat in the form of entropy production as it occurs in the Second Law. 



1.1 Conventional FT 

The first discovery of a FT of the kind we will discuss here was mainly 
numerical in a computer simulation by Evans, Cohen and Morris in 1993 [1]. 
This FT for a non-equilibrium stationary state was inspired by dynamical 
systems theory notions [2] . A similar FT for a transient state was formulated 
in 1994 by Evans and Searles [3]. In this paper we will restrict ourselves to 
the non-equilibrium Stationary State Fluctuation Theorem rather than to the 
Transient Fluctuation Theorem. 

The CFT deals with the fluctuations of work as well as heat in a flnite 
dynamical system. The many particle (Hamiltonian) system is subject to an 
external force, which does work on the system. However, this would heat 
up the system if it were not for an ingenious internal thermostat, realized by 
adding a damping term to the equations of motion. The dynamics of this 
system is purely deterministic. The dissipation taking place in the system 
is manifested in a contraction of the accessible phase space of the system 
which can be related to a (generalized) physical entropy production in the 
system ^ . A FT was derived under a number of assumptions for such a system 
for the heat or entropy production, which will be given below. Some years 
later, Kurchan (in 1998) [6] and Lebowitz and Spohn (in 1999) [7] derived 
a similar FT for a system with purely stochastic dynamics under a number 
of assumptions. 



^ The dynamical systems are required to have an isoencrgetic Gaussian thermostat 
for this to be strictly true, otherwise correction terms appear, though these might 
possibly vanish in the large time or system size limit. 
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In both cases a CFT was found, which can be written in the form: 



P{Qr) 
P{-Qr) 
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(1.1) 



Here P{Qt) is the probabihty that the fluctuating heat produced in the system 
during a time r has a value Qt- Therefore, — Qt is a value of the fluctuating 
heat absorbed by the system of the same magnitude (during an equally long 
time). In taking the limit r — > 00 in eq. (1.1), Qt ~ t scales as r. The 
Qr (and later Wr) in the stationary state can be visualized as fluctuations 
on segments of duration r, obtained by cutting a very long stationary state 
trajectory of the system in phase space into segments. 

For systems with deterministic dynamics the same CFT holds for the work 
done in a time r denoted by Wt, i.e., Qr in eq. (1.1) can be replaced by Wt- 
The reason is that Q^- and Wj- are represented by the same mathematical 
expression here, due to the fact that the thermostat converts all external work 
done on the system into internal heat. 

There is a connection between this CFT and theorems discussed in the 
context of Irreversible Thermodynamics. Gallavotti proved for deterministic 
dynamical systems, that near equilibrium, eq. (1.1) leads to the Onsager rela- 
tions, the Fluctuation Dissipation Theorem and the Green-Kubo relations for 
the transport coefficients [8]. This seems to imply that the fluctuations incor- 
porated in eq. (1.1) go beyond Irreversible Thermodynamics, i.e., beyond the 
linear, near equilibrium regime. There is no estimate available for the range 
of validity of Irreversible Thermodynamics nor are there any results in this 
nonlinear regime to date. 

1.2 The system 

We are interested in this paper in a quasi-many particle system: a Brown- 
ian particle suspended in water and restricted in its motion by a laser induced 
harmonic potential, which is pulled through the water with a constant veloc- 
ity V* (cf. flg. 1). This system was introduced earlier in a somewhat different 
context by Wang et al. in 2002 [9]. 

Contrary to the systems in section 1.1, work and heat are not identical for 
this system. The fluctuations of the work done on it as well as of the heat 
produced by it in a time r were computed based on an overdamped Langevin 
equation [10,11]. Although this system is in principle a many particle system, 
the many degrees of freedom of the water have been contracted to those of 
the (Stokes) friction of the Brownian particle and the strength of its assumed 
Gaussian white noise. In dimensionless units, ^ the Langevin equation reads 



= -it - (xt - X*) + C*. 



(1.2) 



^ Compared to ref. [10], this means k = a = ksT = 1. 
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Fig. 1. Brownian particle in a moving harmonic potential (cf. ref. [9]). 

Here the first term on the right hand side represents the friction of the Brow- 
nian particle in the water, the second term represents the harmonic force due 
to the harmonic potential 



while the third term represents the fluctuations of the Brownian particle due 
to the thermal motion of the water molecules. Xj and x* are the positions of 
the Brownian particle and the minimum of the potential at time t, respectively 
(cf. fig. 1). Essential is that contrary to the pure dynamics of the two previous 
systems, satisfying the CFT, here a mixed dynamics occurs: deterministic due 
to the harmonic force and stochastic due to the water fluctuations. For this 
system, energy conservation reads as the first law of thermodynamics: 



Here Wr is the total work done on the system during time r, i.e., pulling it with 
a constant velocity v* through the fluid over a time r; Q^- is the heat developed 
in the water due to the friction of the Brownian particle during time r and 
AUr = Uri'XT-) — f/o(xo) is the potential energy difference of the particle in 
the harmonic potential in time r. Equation (1.4) clearly shows the difference 
between work and heat in this system. Physically the pulling of the harmonic 
potential drags the Brownian particle along, but with a delay, because of its 
friction with the water, while at the same time this delay necessitates a change 
in its potential energy from its initial position xq to its final position x^. A 
non-equilibrium stationary state will be reached when the friction force cancels 
the attractive force on the particle due to the harmonic potential (cf. fig. 1) 
and the fluctuations around this state will be studied. 

In the non-equilibrium stationary state the averages of the thermodynamic 
quantities and are equal, since AUr = Ut+r — Ut = then. However, 
unlike in the cases of pure (deterministic or stochastic) dynamics, P{Qt) 7^ 



[/,(xO = -(x,-x:)2, 



(1.3) 



Wr = Qr + AUr. 



(1.4) 
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P{Wr), which is the main topic of this paper. 

In the following two sections we first sketch how the distribution functions 
for Wr, Qt and AUr for such a system are obtained, after which the fluctuation 
theorem, involving the ratios of the probabilities of Wt and —Wr, and of Qr 
and —Qt will be stated. 



2 Distribution functions 

2.1 Work 

The probability distribution function of the work WV done on the above 
described system during a time r can be derived directly from the Langevin 
equation [10]. To do that, note first that since the Langevin equation (1.2) 
is linear in x^, the distribution function of is Gaussian. (It is an Ornstein- 
Uhlenbeck process.) Secondly, since Wr — — v* • /J^[xt — Xj] dt depends also 
hnearly on x^, its distribution is Gaussian too. Therefore it is completely 
determined by its first and second moments and can be written in the form: 

P(W^,)~e-(^--'^)'/2.^, (2.1) 

where the variance cr^ = (Wr — WrY- Here the averages are over all time 
segments of the trajectory of duration r. The first and second moment can 
be computed from the solution of the Langevin equation (1.2) [10]. 



2.2 Heat 



A similar simple direct derivation of P{Qt) as was done for P(Wr) cannot 
be performed, since the Qr given by (cf. eq. (1-4)): 



Qr^Wr- AUr 



(2.2) 



is quadratic in x^ via the AUr (cf- eq. (1.3)). A way to obtain nevertheless 
P{Qr) is a much more complicated procedure via its Fourier transform Pr{q). 



PiQr) = ^ J dqPr{q)e-^^'^^. 

— oo 



(2.3) 



Pr{q) can be computed exactly for all r (see ref. [11,13]): 

oo 

Pr{q)= I dQrP{Qr)e"''^^ 



OX]) 



[1 + (1 - e-2^)?2]3/2 



(2.4) 
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Fig. 2. Structure of the complex function Pr{q)e~^'^^^ in the complex plane, with 
p = Qt/Qt = —2, and 4.25, respectively. The saddle points A, B, C and D are 
shown as dots. Wiggly curves along the imaginary axis are branch cuts, ending in 
the branch points qi (top cut) and q2 (bottom cut) indicated with a thin horizontal 
line. The dashed line R is the real axis and Sa, Sb, Sq and Sb are paths of steepest 
descent though the saddle points, at which the function attains a maximum. 

Note that Priq) is a function in the complex plane with branch cuts (due to 
the square root in its denominator) and two singularities at qi = «(! — e"^"^)"^/^ 
and q2 = —Qi, which are also the end points of the branch cuts from qi to 
+ioo and q2 to —ioo (cf. fig. 2). All singularities occur for imaginary g- values 
and introduce exponential (rather than Gaussian) tails in P{Qr) for large Qr- 

The approximate evaluation of P{Qt) can be performed for large r using 
the Saddle Point Method (SPM) [12], capitalizing on the fact that in the 
exponent on the right-hand side of eq. (2.3) is proportional to r. As seen in 
fig. 2, the integral along the real axis R in eq. (2.3) can, for every Q^, 
be deformed to a path of steepest descent S-q in the complex plane that goes 
through the saddle point B, without passing through a singularity — which is 
not possible for the saddle points A, C and D. All that is needed to evaluate 
the integral along the real axis are then the properties of the function at the 
point B. For details we refer to a forthcoming paper [13]. 

To summarize the results of the distribution functions for the fluctuations 
of all three quantities occurring in eq. (1.4): 

1. P{Wr) ~ f.-(Wr-Wrf/2a'^ |g Gaussian, from the Langevin eq. (2.1), 

2. P{/SJJr) ~ e"^'^*^^' is exponential. This can be physically understood by 
observing that for large r, Af/,- is the difference of two independent quantities, 
distributed as the potential energy of a (Brownian) particle in a potential in 
contact with a heat bath, i.e., with a Boltzmann weight e~^^. 

3. Since = — AUr, P{Qt) results from an interplay of the Gaussian 
P{Wr) and the exponential e~^^^'^'^\ Thus the SPM leads to a mixed curve for 
P{Qt) which is Gaussian- like in the center, but has exponential tails, (fig. 3). 
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Fig. 3. Mixed Gaussian and exponential behavior of P{Qr) (sketch). 
3 Fluctuation Theorems 

3.1 Work 

We start by introducing a new formulation of the FT for P{Wr) which is 
more precise than that given (for P{Qr)) in eq. (1.1). Taking the logarithm 
of both sides of eq. (1.1), with Wr instead of Qr, and dividing both sides then 
by Wr, we obtain a reformulation of eq. (1.1) of the form: 

lim^ Fr (pw) = Pw (3.1) 

where Fr is the fluctuation function 



and pw = Wr/Wr is a scaled work fluctuation. Equation (3.1) expresses then 
the CFT for Wr, proven in [10]. A new FT for finite r for the work can be 
obtained from the Langevin equation (1.2) and eq. (2.1), of the form [10]: 



where e(r) = (1 — e^'^)/r. We note that all CFTs used so far in the literature 
restrict themselves to r — oo. In this model one can also discuss the finite 
r behavior and its correction to the infinite time behavior, which shows that 
the slope of Fr{pw) > 1 for all finite r (cf. fig. 4). 

3.2 Heat 

Proceeding similarly with eq. (1.1) for P{Qr) as in 3.1 for P{Wr) one finds 
the CFT in a more satisfactory form: 




(3.2) 




(3.3) 



Jim F,(Pq) =pq 



(3.4) 
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Fig. 4. CFT and NFTs for work and heat fluctuations, 
where pg = Qt/Qt is a scaled heat fluctuation, and 

1 P{Qr) 



(3.5) 



Contrary to eq. (3.1) for W-r, the relation in eq. (3.4) is in fact incorrect, due 
to the interaction of P{Qt) with the (exponential) P{/S.Ur)- 

A NFT can be derived using the SPM [12], both for infinite and for finite r. 
The behavior is determined by the above mentioned singularities in the com- 
plex plane in carrying out the SPM. The result for F^{pq) versus pq is given 
in fig. 4, together with that of Fr{pw) versus pw for comparison. 

While Fr{j)w) versus pw is linear for all r, the behavior of Fri^pq) versus 
Pq is much more complicated. In fact, for r — > oo there are three regimes [11]: 



fApq) 



Pq 

Pq 
2 



for < < 1 
for 1 < < 3 
for Pq > 3 



(3.6) 



Thus, for infinite r, the NFT coincides with the CFT for small fluctuations 
< Pq < 1, then exhibits a parabolic behavior between 1 < pq < 3 and 
finally reaches a plateau where Ft-(pq) = 2 for all pq > 3. The behavior for 
Pq < follows from the asymmetry of F^-^pq) (cf. eq. (3.5)). 

The SPM also allows to study analytically the approach of the finite r 
behavior to the infinite r behavior, giving for large but finite r: 



Fr{PQ) 



Pq + 



2 + 



9{vq) 



+ 



for IpqI < 1 
for Pq > 3 



(3.7) 



where 



8p 



9-p2 



1v* 



In 



(3-p)(l+p) 
(3 + p)(l-p) 



(3. 
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Table I. Results for the work (section 3.1) and heat (section 3.2) fluctuations. 

g{p) = ^8ip-3). (3.9) 

Equations (3.7)-(3.9) show that the asymptotic behavior of Ft{pq) for pq > 3, 
is a slowly increasing function ~ \^Pq^-^, while the asymptotic r — > cxd curve 
is approached as r~^/^ for pg > 3 and as for pg < 1. 



4 Discussion 



In table I, the results for the CFT and the NFT for work and heat are 
summarized. A number of questions and remarks present themselves. 

1. How general is the NFT for systems with mixed deterministic and 
stochastic dynamics? Is the Boltzmann factor and the ensuing exponential 
decay for large fluctuations more general than in this model? One would be 
inclined to think so, in view of the physical argument of point 2 in 2.2. 

2. The relative probability for the Brownian particle to gain rather than 
supply heat to the water is much larger in the NFT than in the CFT (cf. fig. 4). 
This might be of relevance in designing micro or nano machines sensitive in 
their functioning to the heat absorbed during large fluctuations. The CFT 
would not be a good basis to judge this effect. We remark that since Ft{pq) > 
for > in the NFT, > 0, in accordance with the Second Law. 

3. The plateau value of 2 for F{pg) for large pg (> 3) can be understood 
physically. For the probability for large Qr (i.e. pg), the exponential distri- 
bution P(AUr) ~ e"'^'^'^^' dominates over the Gaussian distribution P{Wr) 
(cf. fig. 3). This implies, with Wr = Qr (see the end of section 1.2), that 
P{Qt) ~ e-^l^--^-l so that F^(pq) = 2, with p = 1. 
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4. Similar results as discussed in this paper are obtained for the Transient 
Fluctuation Theorem [3,4]. It is certainly, at least in this model, not the 
identity for all r which obtains in the CFT[10,13]. 

5. All our analytic results have been verified by comparison with the results 
of two numerical methods: a sampling method and a fast inverse Fourier 
transform of eq. (2.3) [13]. In particular, the SPM turns out to give good 
results already for r > 3, whereas curves for r < 3 need to be obtained 
numerically. 

6. The connection between the NFT and theorems of Irreversible Thermo- 
dynamics is unclear, although the linear (CFT-like) behavior of the NFT for 
small fluctuations with < < 1 (cf. fig. 4) suggests that the same relations 
hold as for the CFT for small deviations near equilibrium. 

7. So far, in all cases dealt with here, only one property of the fluctua- 
tions of the thermodynamic quantities work and heat in a non-equilibrium 
stationary state — whose averages are described by the First and Second Law 
of Thermodynamics — has been discussed, viz. the fluctuation function Ft. 
Whether something can be said about other properties of fluctuations of ther- 
modynamic quantities, also beyond the linear regime, remains an interesting 
open question. 
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